Density of states in random lattices with translational invariance 
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We propose a random matrix approach to describe vibrations in disordered systems. The dynam- 
ical matrix M is taken in the form M = AA^ where A is a real random matrix. It guaranties that 
M is a positive definite matrix. This is necessary for mechanical stability of the system. We built 
matrix A on a simple cubic lattice with translational invariance and interaction between nearest 
neighbors. It was found that for a certain type of disorder acoustical phonons cannot propagate 
through the lattice and the density of states g{iij) is not zero at tj = 0. The reason is a breakdown of 
affine assumptions and inapplicability of the macroscopic elasticity theory. Young modulus goes to 
zero in the thermodynamic limit. It reminds of some properties of a granular matter at the jamming 
transition point. Most of the vibrations are delocalized and similar to diffusons introduced by Allen, 
Feldman et al., Phil. Mag. B 79, 1715 (1999). We show how one can gradually return rigidity and 
phonons back to the system increasing the width of the so-called phonon gap (the region where 
g{uj) cx cj^). Above the gap the reduced density of states g(uj)/uj'^ shows a well-defined Boson peak 
which is a typical feature of glasses. Phonons cease to exist above the Boson peak and diS'usons are 
dominating. It is in excellent agreement with recent theoretical and experimental data. 

PACS numbers: 45.70.-n,61.43.Fs,63.50.-x 



In a sharp contrast to electronic properties the nature 
of vibrations in disordered systems is still poorly under- 
stood though these excitations are responsible for such 
important phenomena as specific heat, thermal conduc- 
tivity, propagation of sound and elastic properties. Sohd 
amorphous dielectrics (glasses) are the most common ex- 
ample among these systems. Others are polymers, liq- 
uids, granular media, etc. 

Low frequency plane fong wave acoustical phonons (or 
Goldstone modes) are the most extensively studied delo- 
calized excitations in many disordered materials. They 
propagate through the media ballistically with a speed 
of sound. But we yet do not know what is the upper 
frequency where one can still use this useful notion. One 
can find in the literature hot debates about existence of 
high frequency phonons in glasses. 

Some time ago a new type of delocalized vibrations 
in disordered media was introduced that was called dif- 
fusons [l|. These are vibrations extending through the 
system by means of diffusion. It is an important class 
of excitations which occupy in glasses the dominant part 
of the spectrum. The diffusons may be responsible for 
the thermal conductivity of glasses above the plateau. 
According to 0] the heat in glasses above 20 K is trans- 
mitted by means of a random walk of vibrations from one 
atom to its nearest neighbors. 

In the last years it was discovered that amorphous ma- 
terials as diverse as granular media, foams, emulsions, 
and colloidal suspensions can jam into a rigid, disor- 
dered state where they withstand finite stresses before 
yielding [1, At the point of jamming transition these 
systems are marginally stable but there is a breakdown 
of affine assumptions that underlies their rich mechan- 
ics near jamming. Far away from jamming the density 
of states (DOS) shows Debye-like behavior .g(w) oc w**"^. 



But as the jamming point is approached, both the struc- 
ture of the modes and the DOS exhibit surprising fea- 
tures. In particular g{uj) at low frequencies is strongly 
enhanced and becomes essentially fiat at w = [1]. One 
can conclude that low frequency phonons disappear at 
this point and the macroscopic elasticity theory becomes 
inapplicable. It was found that the low frequency modes 
are neither plane-wave-like nor localized Q and they are 
similar to diffusons 

Recently using a random matrix theory we developed 
a model [8| which we believe can describe typical proper- 
ties of vibrations in disordered systems. Here we present 
an example of a 3d disordered lattice with translational 
invariance where g{uj) is nonzero at a; — > and al- 
most all vibrations are delocalized. We show that affine 
assumptions are violated and no phonons exist in the 
lattice. These properties are similar to those of disor- 
dered systems at the jamming transition point. However 
there is important difference between these two systems. 
We show how one can continuously return rigidity and 
phonons back to the system increasing the width of the 
so-called phonon gap (the region where g(uj) oc w^). 

The vibrational frequencies squared ujf for mechani- 
cal system of N particles are the eigenvalues of the dy- 
namical matrix Mij — / ^rriimj, where $y is a force 
constant matrix and m.^ are the particle masses. The ma- 
trices M and $ are real, symmetric and positive definite 
matrices N x N (scalar model). The last condition is 
important. It ensures mechanical stability of the system. 

To describe vibrations in disordered systems we want 
to use the methods of random matrix theory. One can 
always present every real syinmetric and positive definite 
matrix M in the form 0, [l3| 



M = AA^. 
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Here A is some real matrix of a general form. And vice 
versa for every real matrix A the product AA^ is always 
a positive definite symmetric matrix. In this paper we 
are going to consider disordered systems where are 
eigenvalues of matrix AA^ , with A taken randomly. 
Distribution of eigenvalues of matrix AA'^ was found 
The authors investigated the case where A is 
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a real random matrix with independent elements with 
zero mean (Aij) — and equal dispersions {Af^) — 
(Wishart ensemble ^TS] ) . For a squarte matrix A {N x N) 
the distribution of frequencies g{u!) for N ^ 1 has a 
quarter- circle form (with radius cx ^/N). As a result 
g{Lo) is a constant at w = 0. Matrices AA^ were in- 
vestigated in the theory of financial markets 11 2|. com- 
plex networks 13 1 and wireless communications [ij] . For 
vibrations in random-field spin chains (without transla- 
tional invariance) this approach was used in [lo| . 

In Wishart ensemble each matrix element Mij is in 
general not zero 



ik -^jk • 



(2) 



It corresponds to the infinite-range interaction between 
particles. A similar result for the density of states was 
found for sparse random matrices A with only n ^ N 
non-zero random elements in each row. In the limit N ^ 
I, g{uj) also approaches the quarter-circle form Q 



N 



\/%nr^~^, < w < 2Vy/^. (3) 



but with radius independent of the number of particles 
N. 

Another possibility is to build random matrix A on a 
lattice. Let us consider a simple cubic lattice with N 
particles. Each particle has its unique integer index i 
which takes values from 1 to iV. We build the matrix A 
as follows. The element Aij will have a nonzero random 
value if ith and jth particles are the nearest neighbors 
{i j) or it is the same particle {i = j). All other 
elements Aij — 0. 

In Q it was considered the case of so-called pinned 
latices. Diagonal and non-diagonal elements of matrix 
A were independent random numbers without any corre- 
lation. Thus each particle was randomly connected not 
only with its neighbors but also with space. However the 
Goldstone modes (phonons) cannot propagate through 
such lattice. For existence of phonons it is necessary to 
satisfy also conditions of translation invariance 







(4) 



(we consider rrii — const). It ensures that the Newton 
equations for particle displacements Ui have solutions 
Ui = const for a; = what is the necessary condition 
for existence of Goldstone modes. 
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FIG. 1; The normalized DOS g(i^) in a simple cubic lattice. 
(1) scalar nnm, (2) vector nnm, (3) scalar nnnm, (4) vector 
nnnm (nnm — nearest neighbor model, nnnm — next nearest 
neighbor model (up to the third shell)). Inset: participation 
ratio P{lj) for a) N = 10^ and b) = 27^ 



In this paper we consider a simple cubic lattice with 
translational invariance. We take non-diagonal elements 
of matrix A for nearest neighbors as independent ran- 
dom numbers, but diagonal elements An will satisfy con- 
ditions 



(5) 



Then according to Eq. the Eq. Q will be also satisfied. 

FiglT] shows the normalized DOS g{uj) of cubic lattice 
20 X 20 X 20 with translational invariance (we used pe- 
riodic boundary conditions). The average values of non- 
diagonal elements (Aij) — and dispersion (Afj) — 1 
(for Gauss distribution) . The diagonal elements An were 
calculated using Eq. ([5]). The results are surprising. 
We do not see the expected phonon modes with their 
5p/i(w) oc uj'^ for a; — > 0. On the contrary for scalar model 
with nearest neighbor interactions gicd) even increases at 
small uj. This result is almost identical to g{uj) of similar 
pinned lattice Q. 

The increase of g^td) at small uj is due to weak logarith- 
mic singularity superimposed on a smooth dependence 
<7sm('^)- It exists also for a sparse random matrix A with 
sufficiently small coordination number n p|. However 
with increasing n the singularity is suppressed and g{uj) 
approaches the quarter-circle form, Eq. ([3]), i.e. becomes 
constant at small oj [8j. We verified that flatness of g{u!) 
at small cj is also the case in our cubic lattice if we in- 
crease the number of interacting neighbors or switch from 
scalar to vector model or both (see Fig[T]). Since this 
singularity is not important for the following considera- 
tion we will consider below the scalar model with nearest 
neighbor interaction. 

To understand what kind of vibrations we have in our 
lattice we show on the FiglT] the participation ratios cal- 
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culated according to 
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Here ei{uj) is iih coordinate of the normalized eigenvec- 
tor with frequency u. As can be seen ah modes with 
exception of high frequencies are delocalized. They have 
P{ijo) ~ 0.2 which is independent of the system size. Sim- 
ilar results were obtained for pinned lattices Q. In both 
cases we identified these excitations as diffusons. It was 
shown that the energy transfer through the lattice has 
character of diffusion and the level spacing distribution 
obeys Wigner- Dyson statistics. It also indicates mode de- 
localization. It is remarkable that diffusons turned out to 
be not sensible to the presence or absence of translational 
invariance. These results will be published elsewhere. 




FIG. 2: Non-afSne particle displacements in different atomic 
layers. Solid line is the prediction of the elasticity theory. 



Absence of phonons signals that macroscopic elastic- 
ity theory becomes inapplicable. It is similar to systems 
at the jamming transition point. To check this idea we 
stretched our sample 50 x 50 x 50 for two opposite ends 
with displacements u = ±1 at the ends correspondingly. 
Then we calculated the particle displacements u in each 
atomic layer (2500 dashes). They are shown on Figl2] 
The displacements do not obey the elasticity theory pre- 
dictions shown by the solid line. The fluctuations in dis- 
placements of most of the particles are of the order of 
unity. It means that afline assumptions are violated and 
hence no low frequency phonons (plane waves) can ex- 
ist in these lattices in spite of the translational invari- 
ance. Non affine displacements and violation of elastic- 
ity theory were also found in computer simulations of 
small amorphous bodies formed by weakly polydisperse 
Lennard- Jones beads [lit . 

To measure the Young modulus of the lattice we 
stretched our sample at two opposite ends with forces 
f = ±1 per each particle (at the end) correspondingly. 
The mean sample extension was calculated as difference 
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FIG. 3: The distribution of Young modulus E (multiplied by 
number of particles N) for three different samples. Inset: the 
distribution of elastic constants, (k) = 0.5, (fe^) = 5.5. 



between mean displacements of particles at the ends 



AL = 



Then the Young modulus (in our scalar model) is 



E 



\f\L 
•mil 



(7) 
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where L is the sample length and ao = 1 is the lattice 
constant. First we found that the Young modulus is a 
strongly fluctuating quantity. It depends on the choice 
of the random matrix A. Secondly, the average value 
(E) oc 1/N and decreases with the system size. In the 
thermodynamic limit {N — >■ oo) the average Young mod- 
ulus is zero. It is similar to results of 17| for granular 
matter where static shear and bulk moduli of the disor- 
dered system of interacting particles approach zero at the 
jamming transition point. Fig. |3] shows the distribution 
of Young modulus E multiplied by the number of par- 
ticles N for 3 different samples. All points lay perfectly 
on the straight line (in semi- logarithmic scale) . It means 
that distribution function p{E) has a form 



p{E) = aNe' 



-mEN 



0.36 « 1/e. 



(9) 



The distribution of elastic constants k = —Mi^j is shown 
on the same figure {N = 50'^). There are a lot of negative 
springs (with /c < 0) in the system (about 45%). 

The similarity of our disordered lattice with granular 
media can be also found in distribution of random forces 
in deformed lattice. We have loaded our cubic sample 
20 X 20 X 20 with gravitational forces equal for each par- 
ticle. Then at the bottom of the sample the distribution 
of contact forces was found to be exponential, similar to 
what was measured in a granular matter [18]. 

To get phonons we should introduce a finite rigidity 
into the system. For that let us consider a dynamical 
matrix of the form 



(10) 
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FIG. 4: The normalized DOS g{oj) (solid line) and participa- 
tion ratio P{u)) (dots) for dynamical matrix M = AA''" + fiMo 
with fi — 1. Inset shows the Boson peak in g{Lj)/uj^. 
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FIG. 5: The normalized DOS 3(1.1;) corresponding to dynam- 
ical matrix M = AA^ + j-iMo with different /i. 



Here matrix A is the same random matrix in a 3d cu- 
bic lattice with translational invariance. Matrix Mq is 
a regular dynamical matrix for the same lattice with 
unit masses and all spring constants (between the nearest 
neighbors) equal to unity. 

The results are shown on Fig |3] for = 1 and N = 20^. 
Now g{uj) cx at a; — ^ 0. So we revive phonons at 
small uj below the maximum. The participation ratio 
P{ll!) {N = 27^) in this frequency range has increased by 
more than factor of two (up to the crystal value) what 
is a clear signature of existence of well defined Goldstone 
modes. At higher frequencies (above the maximum) we 
have the same delocalized excitations as on FiglT]i.e. dif- 
fusons. On the inset we can also see that in the transition 
region between phonons and diffusons the reduced den- 
sity of states g{uj)/uj'^ shows a well defined peak. This 
is a well known Boson peak which is a typical feature of 
glasses [1^. We see that in our model phonons cease to 
exist above the Boson peak. This is in agreement with 
results of papers 23-22| derived for glasses. Boson peak 
in our disordered lattice has a non-phonon origin. We do 
not believe it is related to Van Hove singularity shifted 
due to disorder as was found in [23j, [24 1 . 

As a result introducing finite values of /i we open in 
the system a phonon gap — Fig. [S] It is a frequency 
range where g{uj) oc and phonons exist as well de- 
fined excitations. Above the gap phonons cease to exist 
and diffusons are the only delocalized excitations. The 
width of the gap and the frequency of the Boson peak 
increases oc y/Ji. Similar changes of the spectrum take 
place when compressing the marginally jammed solid to 
higher packing fractions [sj. 

Concluding we have shown that for some class of dis- 
order in lattices with translational invariance acoustical 
phonons cannot propagate through the system since the 
macroscopic elasticity theory and affine assumptions be- 
come inapplicable. The major part of vibrations are delo- 
calized diffusons. Their DOS is not zero at w = 0. From 
this point our disordered lattices have similar properties 
with granular systems at the jamming transition point. 



We have shown that this unusual behavior is due to ex- 
istence of high concentration of negative springs which 
makes our lattice extremely soft. It is different from gran- 
ular jammed systems where a coordination number (the 
mean number of contacts per particle) plays an impor- 
tant role. 

Introducing a finite rigidity we return phonons back to 
the system increasing the width of the phonon gap. In- 
side the gap phonons are well defined excitations. Out- 
side the gap the diffusons are dominating and phonons 
cease to exist. In the transition region between phonons 
and diffusons our system shows a well-defined Boson peak 
which is a typical feature for glasses. As a result our 
simple scalar random matrix model is able to reproduce 
typical properties of various disordered systems starting 
from granular matter at the jamming transition point, to 
jammed systems and finally to real glasses. 

We are very grateful to V.L. Gurevich and V.I. Kozub 
for many stimulating discussions and critical reading of 
the manuscript. 
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